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Abstract. We study the high-energy asymptotics of the steady velocity distri- 
butions for model systems of granular media in various regimes. The main results 
obtained are integral estimates of solutions of the hard-sphere Boltzmann equa- 
tions, which imply that the velocity distribution functions f(v) behave in a certain 
sense as Cexp(— r|w| s ) for \v\ large. The values of s, which we call the orders of 
tails, range from s = 1 to s = 2, depending on the model of external forcing. 
The method we use is based on the moment inequalities and careful estimating of 
constants in the integral form of the Povzner-type inequalities. 
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Introduction 

In this paper we address the problem of high-energy asymptotics for solutions of 
kinetic equations used for modeling dilute, rapid flows of granular media. Granular 
systems in such regimes are interesting from a physical point of view, since they 
show a variety of interesting and unexpected properties. They also appear in a 
growing number of industrial applications. Much of the interest to kinetic models 
in this context comes from the fact that such models provide a systematic way of 
derivation of hydrodynamic equations based on the principles of particle dynamics. 
They are also useful for numerical modeling of granular flows. We refer the reader 
to the review papers [TUJ I2H 122] for a general exposition of the subject. 

In dilute flows, the binary collisions are often assumed to be the main mechanism 
of particle interactions. The effect of such collisions is modeled by collision terms 
of the Boltzmann or Enskog type. An important feature of collisions of granular 
particles is their inelastic character: in each collision a certain fraction of the kinetic 
energy is dissipated. This introduces some interesting features in the equations: 
in particular, the only functions on which the collision operator vanishes are delta 
functions corresponding to all particles at rest (or moving with the same velocity). 

To obtain other nontrivial steady states in granular systems, as a general rule, 
a certain mechanism of the energy inflow is required. Experimentally this can be 
achieved, for example, by shaking a vessel with granular particles. In terms of equa- 
tions, several simplified models have been proposed, in the space-homogeneous case, 
which include forcing terms of various types jSU EH] E3- Examples of such terms 
are diffusion (in the velocity space) and Fokker-Planck operators which correspond 
to a physical model of a system of particles in a thermal bath. 

Other important types of problems which lead to similar equations are related to 
self-similar solutions in the homogeneous cooling problem and the problem of shear 
flow [TJHnjEj. m both cases the equations can be transformed, after an appropriate 
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change of variables, to space-homogeneous steady problems for the Boltzmann-type 
equations with force terms that correspond to the negative or anisotropic friction. 

One of the interesting features of granular flows, which has been studied actively 
in the framework of kinetic theory, is the the non-Maxwellian behavior of the steady 
velocity distributions. In fact, experimental data, theoretical predictions and nu- 
merical evidence suggest that typical velocity distributions in rapid granular flows 
have high-energy asymptotics (or "tails") given by the "stretched exponentials" 
exp(— r\v \ s ) with s generally not equal to 2 (the classical, Maxwellian, case), or 
display power-like decay for \v\ large (see [TTH Ho] and references therein). 

The precise form of the asymptotics is determined by several factors, among which 
are the details of the interactions and the forcing models. In the present paper we 
study the model space-homogeneous system with hard-sphere collisions and four 
types of forcing terms. We distinguish between the cases of (i) diffusion (Gaussian 
heat bath), (ii) diffusion with friction (Fokker-Planck type terms), (hi) negative fric- 
tion (obtained in a self-similar transformation in the homogeneous cooling problem, 
and (iv) anisotropic friction which appears in the shear flow transformation. 

We obtain integral estimates for steady solutions using functionals of the form 
(|1.13|) . which indicate that solutions have high-energy "tails" given by the "stretched 
exponentials" exp(— r\v \ s ), with s depending on the forcing terms. We obtain the 
values s = | for the pure diffusion case, s = 2 for the diffusion-friction heat bath, s = 
1 for the negative friction case and s > 1 in the case of the shear flow. Our method is 
based on representing functionals ()1.13|) in terms of symmetric moments, studying 
the infinite system of inequalities satisfied by these moments and using a sharp 
integral form of the so-called Povzner inequalities, similar to the one studied by one 
of the authors [I] in the case of the classical (elastic) space-homogeneous Boltzmann 
equation. We expect that the estimates obtained for the moment inequalities can be 
used for studying the time-dependent moments and the time-evolution of the tails, 
which should be an object of a separate study. 

The problem of high-energy tails for solutions of the inelastic hard-sphere Boltz- 
mann with diffusion model has been studied previously by several authors 17, 
EH! EH Q5| by the methods of formal asymptotic analysis (a formal argument be- 
comes particularly simple if one discards the "gain" term in the kinetic equations). 
The general idea that appears in those papers is that for functions of the type 
h(v) = Cexp(— r\v \ s ), with s < 2, the "gain" term Q + (h,h) in the collision oper- 
ator is a small perturbation of the loss term for \v\ large. Our approach based on 
studying the moments of the collision terms allows us to quantify this idea for the 
solutions of the original problem, without making apriori assumptions about their 
asymptotic behavior. A rigorous analysis of the problem with diffusive forcing has 
been performed in where it was proved, in particular, that steady solutions 
are infinitely different iable and decay faster than any polynomial for \v\ large. A 
lower bound of the type C exp(— r\v | 3//2 ) was also established by using a comparison 
principle. The problem has also been studied numerically by a number of authors 

E3EDI- 

Another series of related results was obtained for the so-called inelastic Maxwell 
models [5J 123 , which are approximate equations obtained by replacing the colli- 
sion rate in the Boltzmann operator by a relative velocity independent mean value. 
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The equations for the Maxwell models simplify significantly by using the Fourier 
transform, and the equations for symmetric moments build a closed infinite recur- 
sive system 012111 . Using the Fourier transform methods, Bobylev and Cercignani 
[6] found solutions to the inelastic Maxwell model with a Gaussian heat bath which 
have high-energy tails exp(— r\v\). In the case of self-similar scaling for the inelas- 
tic Maxwell models, solutions with power-like tails were found [21 El 125] - and it 
was conjectured by Ernst and Brito jTH] that such solutions determine the universal 
long-time asymptotics of the time- dependent solutions in the space-homogeneous 
cooling problem. This conjecture has recently been proved by Bobylev, Cercignani 
and Toscani [7J |H| . 

It is clear, however, that while Maxwell models may give reasonable approxima- 
tions of the macroscopic quantities, the details of the velocity distributions can differ 
significantly from the hard-sphere case. In particular, this is true with respect of 
the high-energy asymptotics which depends crucially on the behavior of collision 
rate for large relative velocities, as can be easily seen from the formal asymptotic 
arguments of the type presented in [THJ H3 EI] ■ 

On the other hand, there is an noticeable gap in the development of rigorous 
mathematical theory, between the Maxwell and hard sphere models. Therefore, 
the aim of this paper is to develop rigorous methods that would allow us to study 
solutions of the hard-sphere Boltzmann equation, with a particular attention to the 
high-energy asymptotics. 

The paper is organized as follows. In Section ^ we present the problem and 
formulate the main results. One of the the most important technical aspects of our 
study is obtaining a precise integral form of the Povzner-type inequalities, which 
we study in Section |21 Section El is devoted to the moment inequalities specific 
to the hard-sphere case. We formulate the inequalities in terms of the normalized 
symmetric moments which appear as the coefficients of power series expansions 
of functional (jl,13|) . We further study the dependence of the inequalities on the 
parameters to find the conditions under which the sequences of the normalized 
moments have geometric growth. Finally, Section 0] presents the proofs of the main 
theorems. 

Most of our inequalities can be used in the time-dependent case, and therefore, 
we begin the next Section by considering the non-stationary Boltzmann equation. 



We study kinetic models for space-homogeneous granular media, in which the 
one-particle distribution function f(v,t), v G M 3 , t > is assumed to satisfy the 
following equation: 



Here Q(f, /) is the inelastic Boltzmann collision operator, expressing the effect of 
binary collisions of particles, and Q(f) is a forcing term. We will consider three 
different examples of forcing. The first one is the pure diffusion thermal bath 
EH HH1 , m which case 
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where \x > is a constant. The second example is the thermal bath with linear 
friction 

&(/) = /iA/ + Adiv(i;/), (1.3) 
where A and \i are positive constants. 

The third example relates to self-similar solutions of equation (jl.ljl for Q(f) = 
I2HIIISI- We denote 

f(v,t) = -±^ f(v(v, t),t(t)), v = ^ y 

where 

1 K 

v (t) = (a + Kt)~ 1 , i(t) = -ln(l + -t), a,n>0. 

k a 

Then, the equation for f(v,t) coincides (after omitting the tildes) with equation 
(jl.ljl . where 

3 (/) = -«div(v/), «>0. (1.4) 
Finally, the last type of forcing is given by the term appearing in the shear flow 
transformation (see, for example, [TT| IH]) 

&(/) = (1-5) 

where k is a positive constant. 

We assume the granular particles to be perfectly smooth hard spheres performing 
inelastic collisions characterized by a single parameter: the coefficient of normal 
restitution < e < 1. To define the collision operator we write 

Q(fJ) = Q + (fJ)-Q~(fJ), (1.6) 

where the "loss" term Q~(f,f) is 

Q-(f,f)=f(f*\v\), (1-7) 
and the "gain" term Q + (f, f) is most easily defined through its weak form: 

Q + (f,f){v)ip(v)dv=±- [ [ f(v)f(w)\u\ ! 4>(v')dadwdv, (1.8) 
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where ii = t> — w is the relative velocity of two particles about to collide, and v' 
is the velocity after the collision. The collision transformation that puts v and w 
into correspondence with the post-collisional velocities v' and w' can be expressed 
as follows: 

v' = v + — (\u\a — u), 

\ (1.9) 
w = w — — (\u\a — u), 

where we set j3 = and < e < 1 is the restitution coefficient. Notice that we 
always have \ < (3 < 1. 

Combining (jl.7j) and (jl.8j) and using the symmetry that allows us to exchange v 
with w in the integrals we obtain the following symmetrized weak form 

Q{fJ){v)^{v)dv = \ f [ f(v)f(w)\u\A P Mv, W )dwdv, (1.10) 

^ ./TCP3 ./lD3 
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where 

AM]{v,w) = ^- I (V>(0 + VK) - 1>(v) - V»H) da. (1.11) 

The weak form ()1.8|) will be sufficient for the purposes of our study. The usual 
strong form |TDJ |2U 123 can be obtained from ()1.8|) by taking ipiy) = 5{v — v ) (see 
also P2|). 

We will assume that the solutions are normalized as follows 

/ f{v,t)dv = l, [ f{v,t)vidv = 0, t = l,2,3. (1.12) 

Since the expected behavior of solutions for \v\ large is C exp(— r\v \ s ), we introduce 
the following functionals: 

FrAf)= I f(v)exp(r\v\ s )dv, (1.13) 
Jm. 3 

and study the values of s and r for which these functionals are finite. This motivates 
the following definition. 

Definition. We say that the function / has an exponential tail of order s > 0, if 
the following supremum 

r* = sup {r > | J>, a (/) < +00} (1.14) 

is positive and finite. 

In the case s = 2 the value of (r*) _1 is known as the tail temperature of / jlj. 
It is easy to see that the number s in the above definition is determined uniquely. 
Indeed, if for certain s > 0, 

< r* < +00, 

then we have r*, = +00, for every s' < s, and also r*, = 0, for every s' > s. 

Another useful representation of the functionals ()1.13|) is obtained by using the 
symmetric moments of the distribution function. Setting 

m p = f{v)\v\ 2p dv, p>0, (1.15) 

and expanding the exponential function in 1)1.13)1 into the Taylor series we obtain 
(formally) 

r / J^L r k \ m sh. 

*■.</>=/„/<•> E ¥ H si )* = E^-*- <" 6 ) 

>/JR fc=0 fc=0 

Then the value r* from p. 14)1 can be interpreted as the radius of convergence of the 
series (J1.16)) . and the order of the tail s is therefore the unique value for which the 
series has a positive and finite radius of convergence. 

We can now formulate the main results of this study. Our first result concerns 
the steady states of equation (jl.l|) corresponding to the first three types of forcing. 

Theorem 1. Let fi(v), i = 1,2,3, be nonnegative steady solutions of the equations 
(jl.ljl . with the forcing terms (|1.2|) . (|1.3J) and (|1.4j) . respectively, and assume that 
fi(v) have finite moments of all orders. Then fi(v) have exponential tails of orders 
|, 2 and 1, respectively. 
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For the shear flow model (jl.5|) . we obtain the following weaker result. 

Theorem 2. Let /4(f) be a nonnegative steady solution of the shear flow model 
()1.5|) that has finite moments of all orders. Then the supremum r\, defined in 
(|1.14J1 . is finite, and therefore, s > 1. 

Remark. The assumption of finiteness of moments of all orders is obviously required 
for the functionals (jl.l3j) to be finite and for the expansions ()1.16|) to make sense. 
However, the moment inequalities we establish below also imply the following apriori 
estimates for all cases of the solutions: Suppose that a moment m Po of any order po > 
1 is finite. Then, in fact, all moments are finite and the solutions have exponential 
tails of the corresponding order. This observation is important, since it excludes the 
possibility of power-like decay for solutions of the considered equations, as soon as 
solutions have finite mass and finite moment of any order higher than kinetic energy. 

The approach that we take in order to establish the above results is based on the 
moment method, in the form developed by one of the authors jl], for the classical 
space-homogeneous Boltzmann equation. We study the moment equations obtained 
by integrating (Jl.lj) against |i>| 2p : 

^=Q P + G P (1.17) 
(in the steady case the time-derivative term drops out), where 

Qv = I Q(f,f)\v\ 2p dv and G p = ! g(f)\v\ 2p dv. (1.18) 

JR3 J R 3 

To investigate the summability of the series (jl.lfij) we look for estimates of the 
sequence of moments (m p ), with p = 4^; k = 0, 1, 2 . . ., and study the dependence 
of the estimates on s. We will be interested in the situation when the series has 
a finite and positive radius of convergence, which means that the sequence of the 
coefficients satisfies 

ITLsk 

cq k <-rj-<CQ\ A; = 0,1,2..., 
for certain constants q and Q > 0. 



2. POVZNER-TYPE INEQUALITIES FOR INELASTIC COLLISIONS 

In this section we establish an important technical result that will allow us to 
control the contribution of the "gain" operator in the moment equations ()1.17|) . We 
consider radially symmetric test functions ip(v) = ^(\v | 2 ). The weak form (jl.8|) of 
the "gain" operator can then be written as 

[ Q+(fJ)(v)y(\v\ 2 )dv = l[ [ f(v)f(w)\u\A+m(v,w)dwdv, 

where 

4M(«,«0 = ^ (*(KI 2 ) + *(K| 2 )) da. (2.1) 

A series of results jSHl ECU I2H C2 EH HI 123 obtained in the case of the classical 
Boltzmann equation develops the general idea that for convex functions \l/ the ex- 
pression (|2.1j) is in a certain sense "smaller" than the corresponding contribution of 
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the "loss" term, which is given by 

A~[V](v,w) = ^(M 2 ) + ^(\w\ 2 ). 

This type of results is generally known as Povzner-type inequalities. An approach 
for obtaining such inequalities in the case of inelastic collisions has recently been 
developed in [TH] . However, for the purposes of the present study we need a better 
control of the constants in the inequalities than those provided by the results of 
[T§] . We will therefore establish a sharper version of the Povzner-type inequality 
for inelastic collisions, using the ideas of [3]. The key point, as in [3], is to look for 
estimates of the integral quantity ([2.1)1 . rather than for pointwise estimates of the 
integrand. 

Lemma 1. For every [3 G [|, 1] there exists a function gpin), on fi G (—1,1), 
which we define explicitly in the course of the proof, such that gpin) is nonnegative, 
continuous, even, nondecreasing for /i G [0,1], satisfies 

2 I g p {2z-l)zdz=l, 
Jo 

and 

^(/i)<l+(-^-l) 2 , 
and for every smooth function ^(x), defined for x > 0, nondecreasing and convex, 

A^)<2 [ gp(2z-l)y(z(\v\ 2 + \v t \ 2 ))dz. 
Jo 

Remarks. 1) The above inequality is a generalization of inequalities (12), (16) 
from [4 to the inelastic case, under the extra assumption of \I/ being nondecreasing. 
Indeed, from the conditions on gp(fJL) it is easy to see that in the elastic case (3 = 1 
we must have gi(fj) = 1. Then the inequality of the lemma reduces to 

A^} < 2 [ V(z(\v\ 2 + \v*\ 2 ))dz, 
Jo 

which is equivalent to the form used in 2) The smoothness assumption on \1/ 
can be dispensed with relatively easily, by elaborating some of the arguments we 
use in the proof. On the other hand, in the most important examples ^>(x) = x p 
with p > 1, which will be used in the moment estimates, the needed smoothness is 
readily available. 

Proof In the proof we use repeatedly the following argument: suppose that ip(v), 
v G M. 3 is a convex function. Then, for almost every a G 1R 3 , and for every b G M 3 , 

ip(a + tb) + if)(a - tb) (2.2) 

is a nondecreasing function of t > 0. (To see this differentiate ()2.2|) in t and notice 
that a convex function has almost everywhere a nondecreasing directional derivative 
in every direction.) 

To apply this argument we notice that since \fr(x) is convex and nondecreasing, 
then also ij){v) = ty(\v\ 2 ) is convex as a function of v G M 3 . In order to introduce 
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the symmetric structure (|2.2|) into the integrand of (j2.1j) we then rewrite v' and w' 

as 

v' = U+-, w' = U--, 
2 ' 2 ' 

where U = is the velocity of the center of mass, and u' is the relative velocity 
after the collision. We then have, according to fll.9|) . 

u' = (1 - (3)u + P\u\a, (2.3) 

where it is the relative velocity before the collision. Further, we pass to the spherical 
coordinates (p, uj) in v! by setting 

u' = p\u\u, p G R, G S 2 . 

Denoting by i/ the unit vector in the direction of it and using ()2.3j) . we can write 

pu = {1 - p)u + /3<t. (2.4) 

We then perform the change of variables from a to w in the integral (j2.1|) . To 
do so, notice that for every test function <p(k), k G M 3 , we can formally extend the 
integration from S 2 to M 3 by writing 

/ ^(a)d<x = f 5( l -^ 7 P±) V (k)dk J (2.5) 

where 5 is the one-dimensional 5-function. 
Changing variables from k to 

k' = (1 -p) u + pk, k' = pu 

and then passing to the spherical coordinates (p, to) we can rewrite the integral ()2.5|) 

as 



i f r ^ - (i -w-ps ,^ - (i -/?),- dpdu 



where 

a= (1-/3)(z/-cj) and 6 2 = 2/3-1. 
The radial integration in ()2.6|) can be performed explicitly, since 



(2.6) 



2 

(a + Va 2 + & 2 ) 2 
\/a 2 + b 2 

After the radial integration, p in ()2.6|) will be expressed through p = (v-u) according 
to 

p = X(p) = (1 - P)p + ^(1-/5)2^ + 2/3-1. (2.7) 
(The last equation is nothing but the condition \k\ 2 = 1 expressed in the new 
variables.) Thus, we obtain the formula 
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where ^ 

9M = p v^+¥ = m^^Wr (2,9) 

Applying identity (|2~51) to (f2~TT) we get 

4s M = ^ ^ 2 ^ ■ + A^f) + - A Y^f) } ^ • ( 2 - 10 ) 

Our next goal is to simplify ()2.10|) to get a convenient upper bound. First, due to 
the convexity of \I/(| • | 2 ), the expression in braces, considered as a function of A > 0, 
is monotonically increasing. Using ()2.7j) it is easy to see that 

< 2/3-1 < A(/i) < 1, 

for all /j, G [—1, 1]. Thus, estimating A by one and setting 

Id 2 

E = 2U 2 + ^- = M 2 + H 2 , (2.11) 

we find: 

f =H!£!M( m . w ) i (2.i3) 



where 



and m is the unit vector in the direction of U. 

We further symmetrize ()2.12|) by using the change of variables uo \— > —a;. Since 
the expression in braces is invariant under this transformation, we can replace the 
function gp(fi) in (|2.12j) by its symmetrized version 

to(ti = \(9M+9p(-IJ))- (2-14) 

It is now a somewhat tedious but straightforward computation to check that gp(fji) 
has the properties listed in the formulation of the lemma. Noticing that 

SH < 1 

E ~ 

and using the convexity argument again, this time for the function ty(E(^)), we 
can replace £ in ()2.12|) by {m ■ uj). 

Next, we see that, for \U\ and \u\ fixed, the integral ()2.12|) has the structure 

<^i(f • uj) ip 2 ( m ■ oj) du , 

where tpi and ip 2 are nonnegative, even and monotonically increasing on [0, 1]. It is 
easy to show that the maximum value of the integral is attained when the vectors 
v and m are parallel. The integral in (J2.12|) is then bounded by 



(2.15) 
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Using that g~p is an even function of \i and performing the change of variables z = 
(1 + /i)/2, we arrive at the conclusion of the lemma. □ 

In the case when the function ty(x) is a power function of x, the bounds of the 
lemma take a more explicit form, and we obtain the following important corollary. 

Corollary 3. Let if){v) = \v\ 2p , where p > 1. Then Ap [if)], given by satisfies 
the inequality 

M\v\ 2p ) < i P (\v\ 2 + \w\ 2 y - \v\ 2p - \w\ 2p , 

where the constant jp, defined by ()2.16|) . is strictly decreasing for p > 1 and satisfies 



Proof. Taking \[?(x) = x p , we can write 



\w 



2p 



A p {\v\ 2p ) < A[[V] - \v\ 2p 
Using Lemma ^ with \fr(x) = x p we get the bound 

^[*]<7 P (M 2 + M 2 ) P , 

where 

7p = 2 I g p {2z-l)z p dz. (2.16) 
Jo 

By Lemma ^ 71 = 1, and so, since z p < z for all z G (0,1), we have 

1p < 1, 

for all p > 1. On the other hand, estimating gp(p) by its maximum 

^(l) = l+(^-l) 2 <2, 



we get 



7 P < 4 / z p dz = . (2.17) 

p+l 







This completes the proof. □ 

Remark. The expression (|2.1fi|) for the constant 7 P simplifies in the cases (3 = 1 
(elastic interactions), when 

2 

p+l 

and in the case (3 — 1/2 ("sticky" particles), when 

p2 p + 1 



1p 



2 p - 2 (p+l)(p + 2) 



In the general case the integrand of (|2.16J) is too complicated to yield an answer 
in closed form. However, the bound (|2.17|) is quite useful for p > 3 and shows the 
correct "inverse first power" decay for p large. 
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3. Moment inequalities 

The estimate of Corollary |3] is a crucial step to obtaining the moment inequalities 
in the form characteristic for the Boltzmann equation with "hard interactions": 
|121 H] . The basic estimate takes a particularly simple form when p is an integer, 
since then the binomial formula can be used to obtain the inequality 

P \\ v \ 2k \w\ 2{p ' k) 



MM 2p ] + (i - 7 P )(M 2p + H 2p ) < % (f) 

k=l ^ ' 



(cf. [4, p. 1189]). In the case of non-integer p a similar result can be obtained, after 
we establish the following estimates of the binomial expansion. 

Lemma 2. Assume that p > 1, and let k p denote the integer part of 21-. Then for 
all x, y > the following inequalities hold 

fco-l / \ 

I r>\ . . . 

,P 



J2 (?) ( x V"* + x p ~ k y k ) <(x + yf -x p -y p 

U — 1 \ / 



(3.1) 



< 

k=l 



Remarks. 1) The binomial coefficients for non- integer p are defined as 
p\ _ p(p-l)...(p- k + 1) k>l _ /p\ | 



kj k\ ~ ' v° 

2) In the case when p is an odd integer the first of the inequalities becomes an 
equality which then coincides with the binomial formula for (x + y) p . 

Proof. The proof will be achieved by induction on n = k p = 1, 2, 3 . . . In the case 
k p = the following inequality is satisfied for — 1 < p < 1: 

(x + y) p - x p - y p < 0. 

Next, for n = 1 and 1 < p < 3, using the above inequality and the identity 

< (x + y) p - x p -y p = / p(p - l)(t + r) p - 2 dr dt, 

Jo Jo 

we obtain 

px py 

(x + y) p -x p -y p < / p(p-l)(t p - 2 + T p - 2 )drdt = p(xy p - 1 + x p ~ 1 y), 
Jo Jo 

which provides the base for the induction. 

Assuming now that the inequalities ()3.1j) are true for In — 1 < p < 2n+ 1, we write 

(x + y) p+2 - x p+2 - y p+2 = / / (p + 2)(p+l)(t + T) p drdt. (3.2) 

Jo Jo 
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By induction hypothesis, the right-hand side of (|3.2j) is bounded from below by 

/ [\p+2)( P +l)(t P + T P )dTdt 

Jo Jo 

+ £ J\p + 2)(p + 1) J2 (*) + tP ~ v ) dr ^ 

and from above by a similar expression with k p — 1 replaced by fc p . Performing the 
integration, using the identity 

+ 2)(p+l) fp\ + 



(k+ 1)0- k + 1) \kj \k + l]' 
and noticing that k v + 1 = A: p +2, we obtain the lower bound for (|3.2j) in the form 



k — 1 

(p + 2)(xy p+2 + x p+2 y) + J2 ( P k + \ 

k=i ^ 



(x fc+1 y p+1 - fc + x p+1 - fc y fe+1 ) 



k=l 

and the upper bound with k p+2 — 1 replaced by k p+2 . This completes the induction 
argument. □ 

We now establish the following bounds for the moments of the collision term Q p 
defined in (11.18)) in terms of moments m p of the distribution function. 

Lemma 3. For every p > 1, 

-m p+ i <Q P <-{1- j p ) m p+ i +%S P , 

where 

^£(IW W +-W (33) 

k=l ^ ' 

and 7 P is the constant from Corollary^ 

Proof. Multiplying the inequality of Corollary El by f{v)f(w) \v — w\ and integrating 
with respect to v and w, we obtain 

Q P <— [ [ f(v)f(w) \v-w\ ((M 2 + \w\ 2 ) p - \v\ 2p - \w\ 2p )dvdw 

2 Jr3 JrS (3.4) 

— (1— 7 P ) / f(v) \v\ 2p / f(w) \v — w\ dvdw. 

Jm 3 

The inner integral in the last term can be estimated as 

f(w) \v — w\ dw > \v\. 



The last inequality follows by Jensen's inequality, since / is normalized to have unit 
mass and zero mean, and the function \v — w\ is convex in w for every v fixed. Thus, 
the last integral term in ()3.4j) is estimated below by 

/0) \v\ 2p+1 dv = m p+1/2 . 
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In the first integral term in (|3.4|) we use the inequality \v — w\ < \v\ + \w\ and the 
upper estimate of Lemma El to get 

\v-w\ ((\v\ 2 + \w\ 2 ) p - \v\ 2p - \w\ 2p ) 



kp / \ 

< 5Z \ ) {\v\ 2(k+1/2) \w\ 2{p ' k) + \v\ 2{p ~ k+1/2) \w\ 2k ). 

k=l ^/ 



(3.5) 



Substituting the estimate ()3.5|) into ()3.4|) and performing the integration we obtain 
the upper bound of the Lemma. 

For the lower bound we use the splitting neglecting the nonnegative Q + 

term and estimating the moments of Q~ in the same way as we did for the second 
integral term in ()3.4|) . This completes the proof. □ 

We next apply the bounds for the moments of the collision terms obtained in 
Lemma El to the steady moment equations 

Q p + G p = , (3.6) 

obtained from (jl.!7|) . Under suitable conditions on smoothness and decay for \v\ 
large of the solutions f(v), the moments of the forcing terms are calculated as follows. 
In the case of pure diffusion (jl.2|) we have 

G p = f(v)/iA\v\ 2p dv = 2fi P {2p+ 1) m H . (3.7) 

In the case of diffusion with friction (J1.3)) . we obtain 

G p = [ f(v)(^A\v\ 2p -XvV\v\ 2p )dv 

= — 2A p m p + 2/i p (2p + 1) m p _! . 

Setting fi = and A = — k in the above identity, we obtain the case of the self-similar 
solutions: (jl.4|) 

G p = 2Kpm p . (3.9) 
Finally, for the shear flow term (|1.5J) . we obtain the inequality 



(3.10) 



Hence, combining the bounds of Lemma H3 with the above identities, we find for 
every p > 1, in the first three cases of forcing the double inequalities 

Gp<m p+l <—^(G p + lp Sp), (3.11) 
2 l — 7 P 

and in the case of the shear flow the one-sided inequality 

m p+ i < — — — {2Kprrip + ~ipSp), (3.12) 

where G p are given by (j3.7|) . ()3.8|) and ()3.9|) . and S p is given by ()3.3|) . 

Notice that since the terms G p and S p depend on the moments of order at 
most p (p — g in the case of S p ), inequalities ()3.11|) and ()3.12|) can be "solved" 
recursively. More precisely, assuming some properties of the moments of lower order 



14 A.V. Bobylev, I.M. Gamba and V. Panferov 

we can use the recursive inequalities to obtain information about the behavior of 
the moments m p , for p large. 

In order to study the summability of the series , it is convenient to formulate 
the moment inequalities in terms of the normalized moments 

where a and b are constants to be determined. Indeed, the coefficients of the power 
series represent a particular case of (|3.13j) . with p = y, a = f and 6 = 1. 

We will therefore study the conditions on a and b under which the sequences of 
normalized moments z p = Zsk have geometric (exponential) bounds. 

We will first look for estimates of the term S p in the moment inequalities (J3.11|) 
and (|3.12|) . expressed in terms of the normalized moments z p . We recall the definition 
of the Beta function 

B(x, y) = jf 1 - *r 1 ds = (3.14) 

which will be used in the proof of next lemma. 

Lemma 4. Let m p = z p T(ap + b) with a > 1 and b > 0. Then for every p > 1, 



where 



(3.16) 



S p < AT(ap + ^ + 2b)Z p , 
Z p = max {z k+1/ 2Z p - k ,z k z p ^ k+ i/ 2 } (3.15) 

l<k<k p 

and A = A(a, b) is a constant independent of p. 

Proof. Substituting ()3.13j) in the expression ()3.3j) for S p we get 

k p , * 

S ? = E [l ( r + \ + b ) ~ k )+ b ) ^+1/2 

k=l ^ ' 

+ T(ak + b) T(a(p - k) + | + b) z k z p - k+1 / 2 
Using ()3.14|) . we can rewrite ()3.16|) as 

k p , s 

T{ap +^ + 2b)J2[ P k )( B ( ak + 7; + b ' a (P- k ) + b ) z k+i/2 Zp- k 
k =i V / 

+ B(ak + b, a{p-k) + ^ + b) z k Zp- k+1/2 ^ ■ 

Next, we estimate the products z k+ i/ 2 z p - k and z k z p - k+ \/ 2 by their maximum Z p , 
obtaining the following bound for the sum in ()3.17j) 



(3.18) 



(3.17) 



Z p (?) ( B ( ak + ^ + b,a(p-k) + b) + B(ak + b,a(p-k) + ^ + b) 
k=i ^ ' 

= Z p jf * sf+^^l - sf- 1 Q {s ak (l - s) a(p - fc) + s a ^- k \l - s) afe } ds. 
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Since the expression in braces depends monotonically on a, we estimate it from 
above by setting a — 1. Further, using the lower bound of Lemma |2l the right-hand 
side of (I3.18J) is bounded above by 

Z p f{s a - +b -\l - sf- 1 (1 - s p - (1 - sf) 

J{) (3.19) 

ds , 



where Xp = if p is an odd integer, and 1, otherwise. Neglecting the negative terms 
in 1 — s p — (1 — s) p and using the definition of the Beta function again, we obtain 
the bound 



Bq + b ,b) 



+ Xp(j y )(B(k p + ~ + b,p-k p + b)+ B(k p + b,p-k p +^ + b)y 



(3.20) 



The first term of ()3.20|) is a constant independent on p; to estimate the second term 
we recall the following asymptotic formula for the Gamma functions pQ : 

lim3^±4^- = l, (3.21) 

p^oo I [p -f s) 

for all r, s > 0. Therefore, taking the first Beta function in the second term of 
()3.20|) for definiteness, we obtain 

p\ „ a <u 1 , r(p + l)r(A; p +f + 6)r(p-A; p + 6) 



r(p + f + 26) T(A; P + 1) T(p -k p + l) 

<Cp 1 -^ b kl +b -\ P -k p ) b - 1 . 



A similar inequality can be obtained for the other Beta function term. It is clear 
now that the second term in ()3.20|) is 0{p~ l ) for p — > oo, and since it also is locally 
bounded for p > 0, it is bounded uniformly in p. Denoting now by A = A(a, b) the 
uniform bound of ()3.20j) we obtain the conclusion of the lemma. □ 

Remark. A more careful analysis of the expression (J3.18j) would allow us to obtain 
a sharper upper bound 

Cp- a Z p (3.22) 
for that expression, at least for 1 < a < 2. Thus, the factor T(ap + f + 6) in the 
estimate of the lemma could be improved to T [ap — | + b) . However, the result 
in the present formulation will be sufficient to obtain the necessary bounds for the 
moments, so we will not pursue the improved estimates based on the bound ()3.22|) . 

We next obtain the simplified inequalities for the normalized moments (j3.13|) . 
Substituting ()3.13|) in the inequalities ()3.11|) and using the estimate of Lemma |U we 
obtain in the case of pure diffusion ()3.7)1 

„ T(ap — a + b) /n „. 2u T(ap — a + b) ,„ _ x 

^ r(a" +f + & ) ^ + 1} ^ * ^ * 1^ T(ap + l + b) ^ + 1} ^ 

lp A r(ap + f + 26) 
+ l- 7p T(ap+l + b) Zp > {6 - 26) 
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for all p > 1. In the case of diffusion with friction (|3.8p . the terms 

2A T{aP + b) vz and 2A T ^ + h ) vz f3 24l 

will be added to the left and the right-hand sides of (j3.23|) . respectively. For the 
shear flow case (|3.10J) we obtain 

^ 2k T(ap + b) lv A r(ap+f + 26) 

l-7,r(ap+f + 6) PZp+ l- 7p r(ap+§ + 6) Zp " ^ 

Using Corollary El for every e > and for all p > 1 + e, the constants involving 
7 P can be estimated as follows: 

1 < z ~ < = K e (3.26) 

1 - lp 1 - ll+e 



and 



7p < i^. (3.27) 



1 - 7p P + 1 
Further, using the identities 

zT(z) = T(z + l) and z (z + 1) T(z) = T(z + 2) 

and estimating 



and 



< c 3 < 2p(2p+l) < ^ 

~~ (°P — a + b) (ap + 1 — a + 6) — ' 

ap + ^ + 26-1 <C 4 P+1 



2 4 ' 

we can reduce the inequalities ()3.23|) to 

T(ap~a + b + 2) ^ „ „ T(ap-a + b + 2) 

°^ T(ap + l + b) ^^^^ r(a P +f + 6) 



T(ap + f + 26-l) 
+ CaKe r(ap+f + 6) Zp ' 



(3.28) 



For the additional terms ()3.24jl appearing in the equation with friction, we use the 
inequalities 

c 5 < — 7T < C 5 
ap + o 

to estimate them as 

r(ap + 6 + l) r(qp + 6+l) 

-Cs^A— - z p and - c 5 A — — — - z p . (3.29) 

T(ap +f + 6) y T(ap + f + 6) 

Finally, for the self-similar solution case we obtain the inequalities 

T(ap + b+l) T (ap + b + l) 



T(ap+1 + b) p ~ £ r(ap+f + 6) p 

T(ap+f + 26-1) 

+ C 4 -ftT £ — — — —TT Z p , 

T(ap +f + 6) 
the last of which is also true in the shear flow case. 



(3.30) 
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We now study the inequalities (j3.28j) - (j3.3U() for the values of a — - corresponding 
to the proposed orders of tails of the solutions. In the case of pure diffusion we take 
a = | and the inequalities ()3.28() take the form 

T(-p + - + 2b - 1) 

c 3 H Zp-i < z p+ i <C 3 K £ /2 z v _ x + C A K £ — 3 3 2 — Z p , (3.31) 
2 1 V3P + 3 + °) 

for p > 1 + e. We notice that if b < 1, the asymptotic formula (|3.21|) allows us to 
control the factor in front of the Z p term in ()3.3H) in the following way: 



r(|p+| + 26-l) 1 

C > K - r(j,' +j + ») % for (X32) 



if we take pi sufficiently large. Inequality (|3.31|) then becomes 

1 
2 

In the case of diffusion with friction the choice a = 1 gives us the inequalities 



c 3 /i Zp-i < z p+ i <C 3 K £ fi z p -! + - Z p , for p>pi. (3.33) 



r(p+ i + 6) 

-C 5 K £ Xzp + c 3 fi z p _! < - i - - 2 p+ 1 < -c 5 A z 



r>-± + 2&) 



+ C3 K e /i + C 4 — — . Z p . 

T(p + b+l) 

(3.34) 

Taking now b < | and choosing pi large enough, we obtain using ()3.21jh 

We can then use (13.31(1 to obtain the following simple inequalities 

C 5 K £ \z p >c 3 [i Zp-! - z p+ i (3.36) 

and 



C5 A z p < C3 /i + - C5 A Z p , (3.37) 



1 
2 

for all p>p\. 

Finally, in the case of self-similar solutions we take a = 2, and (|3.3U() becomes 

T(2p + 2b) 

c 5 kz p < z p+ i < C 5 K £ KZp + Ci K £ T ^ 2p + b+ ^ . (3.38) 
We then take b < 1 and choose p\ large enough to obtain 

Inequality (|3.38(1 then simplifies to 

c 5 kz p < z p+ i < C 5 K £ kz p + \z v , for p>p\. (3.40) 

The second of these inequalities is also satisfied in the shear flow case. 

We have now obtained inequalities for the normalized moments (|3.13(1 in the form 
which is simple enough to be analyzed and which, as we will see below, will allow 
us to prove the results about the tail behavior stated in Section ^ Abstracting 
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now from the precise meaning of the terms in inequalities ()3.33|) . ()3.36|) . ()3.37|) and 
f!3.40|) we can say that they express the balance between the "loss terms" (moments 
of order "gain terms" (terms involving Z p ) , diffusion (moments of order p — 1) 

and friction or the force terms in the shear flow (moments of order p) . We notice also 
that inequalities in the form (|3.28|) . ()3.29|) and (|3.3U|) together with the asymptotic 
formula ()3.21|) can be used to actually derive the values of a for which the series 
(I1.16|) has finite and positive radius of convergence. For the sake of simplicity, since 
these values are already known from the formal arguments, we will not perform 
these computations here. 



4. Proofs of Theorems [T] and [21 

Proof of TheoremUl We will establish the following statement that will imply the 
conclusion of the Theorem (see also the Remark that follows Theorem |2J) . We show 
that for every po > 1 there are positive constants c, q, depending on mo and mi 
only, and C, Q, depending on mo, po and m Po only, such that 

rrisk 

cq k <^-<CQ k , (4.1) 

for all A; > -, where s = | in the case of the pure diffusion, s = 2 in the case of 
diffusion with friction, and s = 1 for the self-similar solutions (jl.4j) . Equivalently, 
we can set a = ~, a = 1 and a = 2 in the respective cases and look for the estimates 

cq p <z p <CQ p , (4.2) 

for all p > 1, with z p defined as in 1)3.13)1 . (The constants in (j4.2j) have to be modified 
to match those in (|4.1j) .) 

Notice that it would be sufficient to prove ()4.2|) for a certain value of b > in the 
definition of z p (|3.13|) . Indeed, since 

c lP ^<^±M<c 2 p^, 

1 (ap + b 2 ) 

changing the value of b in ()3.13|) essentially results in the multiplication of z p by the 
factor Cp bl ~ b2 , which can be compensated for by adjusting the constants in ()4.2|) . 
We fix the value of b < 1 so that inequalities ()3.31|) . (|3.33|) and (|3.34|) are available 
for p sufficiently large. 

The proof of the inequalities ()4.2|) will be accomplished in two steps. The first 
one will be to show that (|4.2j) holds on the initial interval, 1 < p < pi, where p\ 

(dependent on p and b) is chosen so that inequalities (|3.32|) and (|3.35j) hold with 

P _ i-po 

e 2 . 

Step 1: Initial interval. We notice that for 1 < p < pi, the Gamma function is 
bounded both from above and from below: 

< * s r5ir+i) s c °' (4 ' 3) 

where for a > and b > the constants Co and Co depend only on a, b and po- 
Thus, on the initial interval it suffices to estimate m p instead of z p in ()4.2|) . 
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To obtain the desired estimate, we first use Jensen's inequality to derive for every 
< p' < p < p" the inequalities 

{m'/Y <m p < (miry. (4.4) 

Taking p' — 1 and p" = po we obtain the bounds 

cq p <m p <CQ p (4.5) 

for 1 < p < po, with c — C — 1, q — mi and Q = Qo = max{l, m l J Po }. 

Step 1: Pure diffusion. We take e = p ^-, use the bounds ()3.2fi|) and ()3.27j) in the 
moment inequalities ()3.11|) . ()3.7|) and estimate 

Sv < 2 P+1 M D , where M D = max {mjjm„_ u i,m t ,im„j}. (4.6) 

' F F \<k<k p 1 ' + 2 ^2 ' > 

We then obtain, for all p > 1 + e, the inequalities 

2 /I (2p + 1) m p _x < m p+ i_ < 2 K e [Lp (2p + 1) m p _ x + K £ 2 p+1 M p . (4.7) 

Now we see that using (|4.7|) we can extend the bounds (|4.5j) (by augmenting the 
constants q and Q if necessary) to the interval |+£<p<Po + |- Using the 
interpolation inequality ()4.4|) we can then extend the bounds (14. 5 j) to all intermediate 
values po < p < p + \. 

Further, by iterating inequalities ()4.7|) we can cover the interval po < p < pi by 
a fixed number of subintervals of length at most ~, so that finally inequalities (|4.5j) . 
with the constants depending on m , mi, po and m po only, will be extended to the 
whole interval 1 < p < pi . 

Step 1: Diffusion with friction. We argue as in the previous case and obtain using 
(j3.11|) . (j3.8|) the following upper bounds for m p+ i'- 

m p+ i < -2K e Xpm p + 2K ef ip(2p+ l)m p _i + K e 2 p+1 M p , (4.8) 

for all p > l + e. Neglecting the non-positive friction term on the right-hand side 
yields the same upper bounds as in the pure diffusion case. On the other hand, the 
lower bound can be written in the form 

2\pm p + m p+ i>2fip(2p+l)m p _ 1 , (4.9) 

which implies that for every p > 1 one of the following inequalities is true: 

2 \pm p > fxp (2p + 1) m p _! or m p+ i > /ip (2p + 1) m p _i 

Combining the two inequalities and using the interpolation inequality ()4.4|) in the 

second of the cases we obtain 

m p > min | f (p + |) "V-i , (np(2p+ 1) m p _i) 2p+1 j. 

This allows us to extend the lower bound ()4.5j) iteratively to the interval 1 < p < pi. 

Step 1: Self-similar solutions. Using the moment inequalities (j3.11|) . (j3.9|) and 
arguing as above we obtain 

2npm p < m p+ i < 2K £ npm p + K £ 2 p+l M p , (4.10) 

for all p > l + e. Using these bounds we extend ()4.5|) to the interval 1 < p < p\ by 
the same iterative argument as in the pure diffusion case. 
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We now pass to Step 2 and use inequalities ()3.33|) . ()3.36|) and ()3.37|) to extend 
bounds ()4.2|) to all p > 1 by an induction argument. The base of the induction is 
established by virtue of the bounds ()4.5|) and ()4.3|) on the interval 1 < p < p\. We 
further verify the induction step separately in each of the three cases. 

Step 2: Pure diffusion. Our aim is to find the constants q and Q in such a way that 
for every n — 1,2,3 ... , the inequalities ()4.2|) for 1 < p < pi + imply the same 
inequalities for p\ + < p < p\ + |. Thus, assuming ()4.2|) for 1 < p < p\ + 2=i 
we use f!3.33|) to find 

c^cq?- 1 < z p+ , < c(c 3 K pof iQ-l + i) Q^. 

2 2 

Taking q < (03 /i) 3 and Q > (2C3 i^ Po we obtain the inequality 

from which it follows that (j4.2j) is true for pi + r ^ < p < Pi + f • 

S'tep Diffusion with friction. The upper bound case can be treated similarly to 
the previous one, with the difference that inequalities (J3.37|) will allow us to increase 
p by one in each step, instead of one half , as in the pure diffusion case. Thus, for 
every n = 1, 2, 3 . . . , we assume (|4.2|) for 1 < p < p x +n — 1 and obtain using (j3.37jh 

?5 ■ 

We now take Q > to obtain the inequalities 

z p+i ^ C Q p+1 , 

which imply the upper bound (j4.2j) for pi + n — 1 < p < pi + n. 

For the lower bound we see that assuming ()4.2j) to be true for 1 < p < p\ + !2 ^-, 
the inequalities (|3.36|) imply that at least one of the inequalities 

K £ C 5 Xz p > ]- c-iiicqP" 1 or z p+ i > ^c z [icq v ~ x 

is true. By choosing q < min < (- C3 jf) 3 , [ we obtain ()4.2|) for p\ + < 

P < Pi + f ■ 

^tep 2: Self-similar solutions. We use inequalities (13.40)) and argue as in the pure 
diffusion case, assuming for every n = 1, 2, 3 . . . that ()4.2|) holds for 1 < p < P\ + I! ^-. 
We then find 

c 5 ncq p <z p+l2 < (c 5 k e kQ-* + i)cg p+ i 

Therefore, taking q < (05 k) 2 and Q > {2C^K £ k) 2 we obtain (Q for p x + ^ < 
P < Pi + f- 

We can now complete the proof of Theorem by an induction argument. □ 

The above proof contains the proof of Theorem |21 as a special case: indeed the 
inequalities for the normalized moments in the shear flow case coincide with the 
upper inequalities for the case of self-similar solutions. The result of Theorem El 
is weaker than in the latter case, since we were not able to obtain suitable lower 
bounds for the moments in the shear flow problem. 
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Concluding Remarks 

The estimates for the normalized moments that we established certainly deserve 
more attention that we were able to give them in the framework of the steady 
solutions to the kinetic equations. In fact, we hope to return to the problem of the 
time-evolution of the tails by the moment method in a separate paper. Another 
promising direction of study stems from the use of the integral bounds together 
with maximum principles for kinetic equations, in the form suggested by C. Villani 
[32] . There are indications that such methods may yield more precise asymptotics 
(in particular, pointwise upper bounds) for some cases of kinetic equations [THj . 
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